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An asymptotically exact result is obtained for the renormalized phonon energy as a function of 
the on-site Coulomb repulsion Uee in the half-filled Hubbard-Holstein model in the strong-coupling 
region at zero temperature. The result is obtained on the basis of the effective Hamiltonian derived 
in the antiadiabatic region with due attention to the transition from CDW- to SDW-regime with 
the increase of Uee ■ Somewhat contrary to naive expectation, the phonon energy is found to shift in 
the insulating phase and the shift becomes larger in the SDW regime than that in the CDW regime. 



PACS numbers: 71.38.-|-i, 71.10.Fd, 71.28.-|-d, 63.20.-e 



I. INTRODUCTION 



For more than two decades, organic materials have at- 
tracted special attention due to the interesting features 
of a complex system in which charge, spin, and phonon 
degrees of freedom correlate to one another E Similar 
complex features have emerged recently in the ferromag- 
netic perovskite compounds hai-xA^yhiO^ where A rep- 
resents a divalent alkaline-earth element such as Ca or 
Sr.a □ Many people also anticipate that the mechanism 
of high- Tc superconductivity in the copper oxides cannot 
be fully understood without a viewpoint to regard itihese 
materials as a charge-spin-phonon complex system.LrQ 

The Hubbard-Holstein (HH) model provides a useful 
framework to investigate such a complex system. Be- 
sides the kinetic-energy term represented by the nearest- 
neighbor electron hopping integral t, it consists of two 
interaction parts. One is the term to treat the electron 
correlation effects through the short-mnge Coulomb re- 
pulsion Uee as in the Hubbard model.S The other is the 
term to couple an electron with local phonons, leading to 
the phonon- mediated electron-electron attraction —Uph- 
The model without the Uee term is called the Holstein 
model which was proposed to consider the motion of a 
conduction electron coupled with intramolecular vibra- 
tions of ions in a molecular crystal.EI This Holstein model 
has a long history of research and in the last several years 
it has been a subject of intensive studies based Hiai»ly on 
the quantum Monte Carlo (QMC) simulationsJljilj 

In the adiabatic region in which t is much larger than 
the local-phonon energy ujq, the HH model as well as the 
Holstein model can be studied vdth-use of the traditional 
Migdal-Eliashberg (ME) theoryJlil'y In the antiadiabatic 
region {t ^ loq), however, a more careful treatment is nec- 
essary and numerical approaches such as an exajc±_diag- 
onalization methodO"E2l and QMC simulationsEilEj have 
been employed. These studies have revealed that the HH 
model can afford interesting phase diagrams represent- 
ing the competition among charge density wave (CDW) , 
spin density wave (SDW), and superconductivity as we 



change the parameters involved in the system, i.e., Uee, 
Uph, t, loq, and n the electron filling. Crudely speaking, 
the CDW (SDW) nature manifests itself in the ground 
state at half filling, if the effective electron-electron inter- 
action U defined hy U = Uee — '^Uph is negative (positive). 

In this paper, we focus on the renormalized local- 
phonon energy ujo in the HH model in the antiadiabatic 
region at zero temperature (T — 0) and inquire how much 
change can be seen in cjo at the CDW-SDW transition. 
Note that the shift Sljq defined by Sloq = loq — Cjq is di- 
rectly connected with density fluctuations the amount of 
which depends sensitively on the nature of the ground 
state. We shall give an exact result for Sloq as a function 
of Uee in the strong-coupling region {Uph ^ wo) and dis- 
cuss its change at the CDW-SDW regime boundary. We 
consider the half-filling situation at which we presume 
that the largest shift in ujq appears. 

The antiadiabatic condition is usually fulfilled in the 
strong-coupling region due to the small-polaronic band- 
narrowing effect. This-iSpalso emphasized by Alexandrov 
and Capellmann (AC)c3'E3 who, starting with a comment 
on the importance of antiadiabatic motion of electrons in 
a strongly-coupled electron-phonon system, have derived 
an analytical result for 5loo presumably at C/ = by 
using the Green's-function method. The result of AC is 
obtained at T larger than the critical temperature for the 
ordered state of CDW or SDW. In such a temperature 
range, the amount of density fluctuations differs much 
from that at T = 0. Thus we expect to obtain a result 
for 5ujQ which is different from that of AC. 

Subtlety exists in the present problem as for the in- 
clusion of an ordering nature of electrons in the Green's- 
function method. In order to avoid the difficulty origi- 
nating from this subtlety, we employ an approach based 
on the effective Hamiltonian 7?eff derived from the HH 
model in the antiadiabatic and strong-coupling region at 
half filling. We observe that H^a for C/ < (the CDW 
regime) is nothing but an Ising-like pseudospin Hamil- 
tonian with an exchange interaction depending on both 
\U\ and local phonon number. On the other hand, Hcs 
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for U > (the SDW regime) can be mapped onto the 
Heisenberg spin Hamiltonian with the same exchange in- 
teraction. 

By evaluating ljq with use of these i?eff 's, we find that 
dujQ depends on the sign of U. If we compare them at the 
same value of \U\, Sujq for [/ > is larger than that for 
U < 0. This result may be understood if we remember 
the fact that phonons with energy much larger than the 
characteristic energy of small polarons or bipolarons can 
react to both spin and charge fluctuations in the elec- 
tronic system. In the HH model in the strong-coupling 
region, spin degrees of freedom give larger fluctuations 
than those of charge. 

In Sec. H we give the Hamiltonian for the HH model 
and explain an overall feature of Qq as we change the pa- 
rameters in the system in order to specify the quantity to 
which we pay attention in this paper. In Sec. Ill , we show 
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various effective Hamiltonians, Hcs's, depending on the 
nature of the ground state in the antiadiabatic region. In 
Sec. IV, the phonon energy shift Suiq is calculated at half 
filling based on the -ffcff's. Finally in Sec. we summa- 
rize our results and briefly discuss a remaining problem 
such as the dependence on n. Supplementary discussions 
are given in Appendices A, B, and C. Throughout this 
paper, we employ units in which ?i = fee = 1. 



II. HUBBARD-HOLSTEIN MODEL 

Consider a system of electrons interacting with one 
another through the short-range Coulomb repulsion and 
hopping around in a crystal composed of molecular units, 
each of which provides an intramolecular high-energy op- 
tic phonon coupled to the electrons. Assuming that both 
the energy of phonon loq and the electron-phonon cou- 
pling constant are independent of wave vectors, we can 
write the Hamiltonian for the A^-site system as 



i,j.{T i i 

+ VauJo pijat + 4) + cjQ alaj, (2.1) 

i i 

where Nij is defined as unity for a nearest-neighbor-site 
pair, and zero otherwise, Ci^ the operator to annihi- 

late a spin-tj electron at site i, riio- = cl^Cic^, pi = ni|+ni|, 
p the chemical potential, a the nondimcnsional electron- 
phonon coupling constant defined through a = Uph/i^a, 
and tti the phonon annihilation operator at site i. 

Information about the phonon properties is involved in 
the phonon Green's function £'(q, ia;^), defined by 

D{q,iLO,) ^ - J^^ dTe^™^(T,0q(r)0_q), (2.2) 

where oj^ = 2'kvT is the boson Matsubara frequency 
with an integer v, Tt is the time ordering operator. 



(j>q = aq + alq with aq = N J2 
the ion-displacement operator, and (• • •) denotes the ther- 
modynamic average. The renormalized phonon energy 
djo appears as a pole of the retarded form of D, D^{q, uj), 
in the complex uj plane. 

If we employ the Green's- function method for the eval- 
uation of -D(q, iuj,^) in the antiadiabatic {t <^ wq) and 
strong-coupling (Uph ^ ujq) region, we encounter a re- 
ally complicated situation because we A^ed to make a 
proper inclusion of the vertex corrections.tj On the other 
hand, an exact diagonalization method for the study of 
a small-size HH system is useful in studying an overall 
qualitative feature of c^o in this region. The size effect 
is expected to be small due to the infrequent hopping 
nature of electrons.E3 In Fig. |l|, we show a typical result 
for the excitation energies in the energy region less than 
3wo as a function of Uee with a fixed value of Uph in a 
two-site half-filled HH model. We obtairL_Lhis result by 
using an exact diagonalization methodllj'Ej to calculate 
the pole positions in D^{q,uj) at q — w (the antibonding 
molecular orbital state) with the parameters t — 0.2luo 
and Uph — 2luo. 

There are basically two kinds of ion displacements in 
this region. One is the shift of equilibrium ion posi- 
tion associated with the change in local electron num- 
ber. The other is the vibration of small amplitude around 
each equilibrium position. Accordingly, ujq has two main 
branches of modes in Fig. 0(a) except for [/ « 0. The 
former mode of purely electronic origin repr|esented by 
the dashed curve was discussed by Ranningeir3 as a lat- 
tice deformation wave driven by charge fluctuation of the 
system. Thus we denote the mode as "charge". (For 
comparison, we have also shown the spin-excitation en- 
ergy by the thin dashed curve in the figure, though it 
does not give a pole structure in D^.) The latter mode 
is an intrinsic phonon mode always appearing at the en- 
ergy slightly less than ujq. We plot this mode by the solid 
curve and label it as "phonon" . 

The charge mode has a large spectral weight of the 
order of 8a in the region of J7 < (the CDW regime) 
as shown in Fig. [l](b) and this is the mode which evolves 
from the softened mode in the adiabatic regionEZl with the 
increase of Uph at a fixed value of Uee- However, it does 
not deserve further investigation, because in the CDW 
regime (Do is known to be vanishingly small, i.e., of the 
order of the exponentially small bipolaron hopping en- 
ergy and in the SDW regime {U > 0) its spectral weight 
is negligibly small. 

The intrinsic phonon mode has the spectral weight of 
about unity irrespective of U. In the following, we mainly 
concern ourselves with the amount of the energy shift Sujq 
of this mode. Note that this phonon mode couples with 
the charge mode at Uee ~ "^Uph + ujQ, originating from the 
coincidence of the bare phonon energy with the charge 
excitation energy at this value of Uee- We will leave a 
detailed study of this special situation in the future. 

Three comments are in order concerning the results on 
the small-size HH model. (1) For the special situation of 
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[/ « with not so strong Uph in which a metaUic state ap- 
pears as a jGonsequence of the competition between CDW 
and SDW^J the third mode plotted by the dotted curve 
in Fig. has also a considerable spectral weight. This 
mode represents the coupling between charge and intrin- 
sic phonon modes. For this reason, we denote this mode 
as "chargc-phoncjail in Fig. ^ (2) The Green's-function 
approach of ACatj involves only two-site processes. The 
lattice effect is included merely through the coordination 
number z. Thus their result should be consistent with 
our two-site calculation on Stoo. We find that the pa- 
rameter dependence, i.e., Sluq oc zt'^ / {o^uiq), agrees, but 
the proportionality constant is different. We also note 
that according to our rigorous argument in Sec. IV, Sloq 



depends on the lattice structure in a more complicated 
way than merely through z even in the strong-coupling 
limit. This is especially the case for the SDW ground 
state. (3) A rather extensive QMC calculation was done 
by Berger et al. on the 4x4 two-dimensional (2D) HH 
model. Their result on Cjq as a function of is given in 
Fig. 8 in Ref. ^ for the half-filled system at q = (p-nTr) 
and with the values of Uph = t = ojq and T = O.lcjo.cS It 
seems that they have detected only the dominant mode 
in each region, i.e., the charge mode for [/ < and the 
phonon mode for U > 0. With such an interpretation 
of their QMC result, we find that they have obtained a 
very similar result to that in Fig. This is encourag- 
ing, because the theory for uiq in the strong-coupling and 
antiadiabatic region seems to be applicable even to the 
medium-coupling, i.e., Uph ~ wq, and not very antiadia- 
batic, i.e., t w wq, case. 



III. EFFECTIVE HAMILTONIAN IN THE 
ANTIADIABATIC REGION 

A. Lang-Firsov transformation 

The Lang-Firsov (LF) canonical transformationil is 
useful to treat the HH model in the strong-coupling re- 
gion. It converts H into H as 



H ^ e^H e-^ = T + V, 



(3.1) 



where P is defined as P = — — Pi(al — T 

is the inter-site hopping term, described as 

f = ~tJ2N.Aa^,.Xjx„ (3.2) 

with Xi = exp[y/a{ai — a|)], and V the on-site term is 
given by 

V = U^ Ui^riii + Wo X! ^ (m + Uph) ^ Pi- (3.3) 

i i i 

Note that the on-site instantaneous interaction is 
changed from Uee into U — Uee~'2Uph- Note also that Cia 



in H should be interpreted as an annihilation operator of 
a "polaron", i.e., an electron dressed with ion displace- 
ments, rather than a bare electron in H. Similarly, 
in H represents an annihilation operator for the quan- 
tized ion vibration around \fapi the shifted equilibrium 
position depending on local electron number. In accor- 
dance with this change of a^, /^(q, zu;^) in Eq. (2.2) can 
be decomposed into three parts as 

D(q, iw^) = Dphisi, iuJu) - 2y/aDcphiq, i^^v) 

-I- 4aDc(q,it^^). (3.4) 

Here the intrinsic phonon part Z?p/i(q, ioju) is the Fourier 
transform of — (Pr</'q(''")0-q)5 with the subscript H rep- 
resenting that the r dependence of 0q as well as the 
thermodynamic average should be considered with re- 
spect to H. The charge part _Dc(q, zwi^) and the charge- 
phonon coupling part Dcph{^^,i^l?) are, respectively, ob- 
tained by the Fourier transform of — (T^- Pq{T) p-q) fi and 
-(Pr('?!'q(T)p-q + Pq (T)(/)_q)) ^ , wlicrc thc dcusity Oper- 
ator pq is defined as pq = N^^/^ J2i e^'^'^'pi. 

In terms of a complete set of energy eigenstates 
satisfying H\m) ^ |m) , we can give the spectral repre- 
sentation for these quantities. For example, I?p/i(q, iWj/) 
is expressed as 

Dph[(\,iuJu)= I -dnj—^^—^Bph{<l,^l.), (3.5) 



TT [iuj^y 



where the phonon spectral function Bph is defined as 



Bph{q,n)=7r{l-e-''^^) 



1 



Err 



-Err,/T 



\{m\<j>q 



^6{Q - Err + E„ 



(3.6) 



At T = 0, we need to cons ider only the ground state |0) 
in thc m sum in Eq. (3.6) and ljq can be calculated by 
the difference between Eq and £"„ for excite states, |n)'s, 
satisfying (O|0q|n) ^ 0. Thus our problem is reduced to 
the analysis of energy eigenstates. 



B. Energy spectrum structure and second unitary 
transformation 

For a further analysis of the HH model in the antiadi- 
abatic and strong-coupling region, we need to know the 
characteristic energy structure of the model. The struc- 
ture can be best explained by starting with the atomic 
limit, i.e., t = so that H — V, and the case of J7 > 0. 
Irrespective of the electron filling n, the total energy of 
an energy eigenstate can be specified by a couple of non- 
negative integers Nph and Nd which represent, respec- 
tively, the total numbers of phonons and doubly occupied 
sites in the state. Of course, a high degree of degener- 
acy exists for the energy ujoNph + UNd— {n + Uph)nN 
due to the many possibilities of choosing sites to place 
up- or down-spin electrons and those to excite phonons. 
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However, the ground-state manifold characterized by 
{Nph, Nd) = (0, 0) is well separated from the first-excited- 
state manifold by the energy of either loq or C7, whichever 
smaller. With the increase of t, a "band" will be formed 
to lift the degeneracy to some extent in each manifold, 
but as long as t is small compared to U or luq, the pair 
{Nph,Nd) — (0,0) remains to be a set of good quan- 
tum numbers for the ground-state "band" . The same is 
true for other excited-state manifolds, if there is no ac- 
cidental degeneracy among different sets of {Nph^NdYs. 
(The degeneracy may occur if the ratio U /uiq is some 
rational number.) Thus, in general, the present system 
has a structure for the energy spectrum such as shown 
schematically in Fig. |2|(a). 

A similar consideration can be made for the case of 
U < 0. Here the ground-state manifold has the largest 
number of under the given total electron number nN. 
Then, if we reinterpret Nd in the above argument as the 
reduction in the number of doubly-occupied sites from 
the ground-state manifold and take \U\ instead of U, we 
obtain exactly the same energy spectrum structure as 
shown in Fig. ^(a). 

A special treatment is necessary for the case of t7 = 0, 
or \U\ ^ < if we take the "bandwidth" into account. 
This is one of the cases in which the ratio U/ujo is equal 
to some rational number. Actually, this has a simpler 
energy spectrum structure than before as illustrated in 
Fig. g(b) , because Nd is no longer meaningful and Nph is 
the only good quantum number here. 

The energy eigenstates in each manifold can be de- 
scribed by an effective Hamiltonian Hcs which can be 
derived from H by eliminating the polaron hopping ef- 
fects on inter-manifold processes, i.e., the virtual transi- 
tions accompanied with the change in the energy iOQNph + 
\U\Nd. For the ground-state manifold, a considerable 
number of literature have already been accumulated, 
to the derivation of for the Hubbard papdel,t2rE3 
the Holstein model,EJO and the HH modeLa Roughly 
speaking, all these theories make a degenerate perturba- 
tion expansion with respect to t up to second or fourth 
order. 

We extend such an approach even to the case of 
excited-state manifolds. For the purpose, we transform 
H into Hcs by introducing some adequate unitary oper- 
ator e^. We define the antihermitian operator R pertur- 
batively on the principle that we retain only such terms 
as to make the energy uj^Nph + \U\Nd invariant among 
the polaron hopping processes in each order of t. 

Following the principle faithfully, we obtain an expan- 
sion form of i?eff up to second order in t as 



cff 



V + h['J + H, 



r(2) 

cff ' 



(3.7) 



where describes the first-order polaron hopping pro- 

cesses and -ff^g represents the polaron-polaron interac- 
tion part induced by second-order virtual hopping pro- 
cesses. In Fig. H, all the processes contributing to 



are shown schematically. Mathematical details to obtain 
and are given in Appendix A. 



C. Spin versus polaron liquid 

Accorc^g to the Holstein's argument on small 
polarons,Ll we readily find that the characteristic energy 
for -ffg^g , ei, is the polaron hopping energy te~". Some 
analysis is needed to find the energy scale to character- 



(2) 

ize H^ff \ 62- The matrix elements related to the hop- 
ping processes among adjacent three sites [Fig. ||(c)-(e)] 
are exponentialhi, small in the antiadiabatic and strong- 
coupling region.oln addition, the matrix element for the 
process (b), the bipolaron hopping process between near- 
est neighbor sites, is also found to vanish exponentially. 
Thus only the process (a) survives. The estimation of-the 
corresponding matrix element suggests £2 = i^/C/p/i.Ej 

Comparison of ei with £2 implies two situations: 
One is the case of ei ^ £2, or equivalently 
{Uph/ujo)e~^^'''^/'^°^ (i/wo), in which we can neglect 



(1) 



("2") 

and need to treat only H^f^ ^ leading to some spin- 



model Hamiltonian for Heg. Thus we call this case "the 
spin region". The other is the case of ei ^ £2 in which 

contributes much to the properties of the system. 
We call the latter situation "the polaron-liquid region" . 

Although the condition to reach the polaron-liquid re- 
gion is compatible with the condition of antiadiabaticity 
and strong coupling [{t/ujo) <C 1 ^ (t/p/j/cjo)], the spin 
region is realized eventually for extremely large values 
of a = Uph/ujQ. Thus we mainly consider the spin re- 
gion in relation to the phonon properties in this paper. 
However, the polaron-liquid region, especially its ground 
state, is quite intriguing in connection with superconduc- 
tivity. Therefore we also give a concrete form of iJoff for 
the polaron-liquid ground-state manifold in Appendix B 
in order to provide a theoretical basis for a future study 
of superconductivity. 



D. Spin Hamiltonian 

We now provide a convenient form of H^s to determine 
uiQ for the intrinsic phonon mode in the spin regime. Note 
that loq is given through the evaluation of Bpfi{(l,i^) in 
Eq. (3.6). The evaluation can be done using Hcs for 
the manifolds labeled by Nph = 0, 1, 2, • • • and Nd — 0. 
Thus, although no restriction is imposed on Nph, we con- 
fine ourselves to the manifolds specified by Nd — 0. As 
for n, we consider only the half-filling case for simplicity. 
Since a concrete form for HaSjVph effective Hamilto- 
nian for the manifold labeled by {Nph,Nd = 0) depends 
critically on the nature of the ground-state manifold, we 
need to treat the three cases, i.e., U > 0, U < 0, and 
U = 0, separately. In the following, we give our final 
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result for HcS-.N^h each case. Some tips for deriving it 
are explained in Appendix A. 

Let us start with the case of C/ > (the SDW region) , 
in which each site is occupied by one electron in the man- 
ifolds with Nd — 0. Suppressing the chemical-potential 
term, we obtain 

^Nph0^o + Y,HhJ ■■ SDW) 

(hi) 



EE 

p,q L 



i-2a)P (-2a)9 
p\{p + Ly.ql(q + L)l 



(3.8) 



where the symbol ^ indicates the sum over nonnega- 
tive integers p and q under the condition of < p -f- g < 
Nph for a given number of Nph and another symbol 
specifies the sum over integer L from — min(p, q) 
to infinity provided that the term L = — |?7|/ljo should 
be skipped, if |?7|/a;o is an integer. The spin operator 
h{i,j : SDW) is defined as 



/i(z,j:SDW)-2(S, -S, -p,p,/4). 



(3.9) 



with Si • S, 



- Sf ^1 + {S+S- + 5r^+)/2, where = 

(n^i -n.ii)/2, 5+ = c|^c^, and 5*^^ = Although 
the presence of phonons brings about some modification, 
it seems to be quite natural to obtain the Hcisenberg- 
like spin Hamiltonian in this SDW region in which the 
electron spin resides at each site. 

Difference appears in the case of C/ < (the CDW re- 
gion), in which the ground-state manifold has the largest 
number of doubly-occupied. In this case, we obtain 
HcS:Nph in the form of Eq. (^^) with the replacement 
of the spin operator h{i,j : SDW) by h{i,j : CDW), 
defined by 



1/4) 



(3.10) 



where rjf r-iiPi ~ l)/2 is the z-component of pseudospin 
operators.ES Note that the pseudospin transverse com- 
ponents corresponding to the bipolaron hopping process 
[the process (b) in Fig. |j vanish because of the immobile 
nature of bipolarons. 



At [/ = 0, H, 



is also obtained in the form of 



Eq. (^), but h{i,j : SDW) should be replaced with 
h{i,j : U — 0), defined by 

h{ij : [/ = 0) = h{i,j : SDW) + 2h{iJ : CDW) 
+ (p. + Pj)/2 

= 2(S, • S, + p,p,/4) - (p, + p,)/2. (3.11) 

This expression visualizes perfect symmetry of the sys- 
tem with U = between spin and charge fluctuations. 
This result at Nph = was derived in Ref. |3^ . 

With these JJ-dependent spin operators h{i,j : U), 
-ffcffio the Hamiltonian to describe the ground-state man- 
ifold is reduced further as 



(3.12) 



where Jq is defined as 



Jo(f/ee,C/..) = 2t^e-f:i^^^, (3.13) 



for U ^ 0. At f/ = 0, the £ sum starts with £ = 1. For 
large a — Uph/ojQ, we can evaluate the ^ sum easily and 
obtain 



2^ 



\U,e - 2Uph\ + 2Uph 



(3.14) 



For J7 > or Uee > 2Uph, this is reduced to 2t'^/Uee, the 
Anderson's super-exchange interaction in the Hubbard 
model.l23 



IV. PHONON ENERGY SHIFT 

A. Ground state 

With reference to the conventional theories on the spin 
Hamiltonians, we readily know the ground state |0) of 
HcS;o- For U > 0, the antiferromagnetic (AF) state in the 
Heisenberg model appears as the ground state with Eq 
the ground-state energy depending on the lattice struc- 
ture. In general, Eq per site is given as 



Eo/N = -zJoS {S + S)~ z Jo/4, 



(4.1) 



where S{— 1/2) is the spin magnitude and 6 represents 
the amount of spin fluctuations in the AF state. In the 
two-site case with z = 1, S is equal to unity. Another 
exact result is known for the one-dimensional (lD)|-chain 
with z = 2, in which S is equal to 21n2 - 1 = 0.386EO For 
higher dimensions, Eq is not known exactly, hut the spin- 
wave theory provides reliable estimates for For 2D 
square [z = 4], simple cubic [z = 6], and body-centered 
cubic (bcc) [z — 8] lattices, the values of S have been 
estimated as 0.158, 0.097, and 0.073, respectively. 

For U < 0, the "AF" state also appears as the ground 
state, but "spin" fluctuations are absent due to the Ising 
nature. Thus we obtain Eq in the form of Eq. ( [4.l[ ) with 
S = irrespective of the lattice structure. Physically, the 
present "AF" state represents the bipolaron-ordered bi- 
partite lattice (the CDW state) and the absence of "spin" 
fluctuations reflects immobility of the bipolarons. 

At [/ = 0, we should consider the competition of these 
two states. Although the energy difference is small, it 
turns out that the Heisenberg-type AF state appears as 
the ground state with the same amount of spin fluctua- 
tions as that for U > 0. This is because the transverse 
spin components in the effective Hamiltonian for U — 
are just equal to those for U > 0. Thus Eq is given in 
the form of Eq. (|1|) with Jq at Uee = '2Uph in Eq. ( |3.13| ) 
and the same d depending on the lattice structure. 
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B. One-phonon excited states: Numerical study 

The effective Hamiltonian for the manifold specified by 
Nph = 1 and Nd = is given by 



-ffeffii = Wo + ^ h{i,j : U) 

X [Jo + Ji(a| - a]){ai ~ aj)/2], 



(4.2) 



with Ji, defined as 



2a 



\U\-Lua 

{2a-i-\f _ 
^+1 



E 



1 



|C/|+^c^o 



(4.3) 



for [/ 7^ 0. At [/ = 0, Ji should be evaluated with the i 
sum starting with I = \ instead of ^ = 0. In the strong- 
coupling limit, J\ is reduced to 



Jl{Uee, Uph) 



2U, 



ph 



aUph \ \Uee - 2Uph \ + 2Uph 



(4.4) 



Comparison between Eqs. ( ^.14 ) and (4.4) indicates that 
Ji/Jo is small, i.e., of the order of 1/a. 

Now we investigate the intrinsic phonon mode which 
appears as a peak in the spectral function i?p/i(q, 51) in 
Eq. ( [3.6| ). The one-phonon excited state \n) giving rise 
to the peak should not only be an eigenstate of -ffoffii but 
also satisfy the property of 



1. 



(4.5) 



In pursuit of such an eigenstate, we have examined all 
the eigenstates of iJcffii for small-size ID chains with the 
periodic boundary condition by the exact diagonaliza- 
tion method. This examination clarifies th at there ex- 
ists one eigenstate \'n{q)) satisfying Eq. (4.5) for each q. 



(Other eigenstates give spectral amplitudes of the order 
of at most t'^.) Using the corresponding energy eigen- 
value En{q), the energy shift for the intrinsic phonon is 
given by Sujo{q) = ujo + Eq - En{q). 

In Fig. ^, an example of the calculated results for 
duJo{q) in units of Jo is given for the ID chain up to 
ten sites. We observe some characteristic features for 
Sujo{q) from this calculation: (l)For U < 0, the size de- 
pendence for SujQ{q) is negligible, suggesting that this 
quantity is determined only by the processes between 
nearest-neighbor sites. (2)For U > 0, the size dependence 
is appreciable, indicating the importance of the longer- 
range fluctuation effects. (3)The SDW regime {U > 0) 
has a larger 6ujo{q) than the CDW regime {U < 0) for 
any value of q, implying that the former provides larger 
density fluctuations than the latter. A physical inter- 
pretation of this result is the following: The amount of 
spin fluctuations is larger than that of charge fluctuations 



in the strong-coupling region. Now, for the high-energy 
phonons, even spin fluctuations can be detected as a kind 
of density fluctuations of electrons. Thus the SDW phase 
gives larger 6uJo{q)^s than the CDW phase. (4)Irrespec- 
tive of the sign of U, the function form of Su!o{q) is seen 
to be 1 — cosq. 



C. Asymptotically exact expression 



A closer examination of the above numerical study of 
i?eff:i reveals that we can write \n{q)) as 



l"(q)) =<|o), 



(4.6) 



with |0) the ground state of J?cff:0- Although this is an 
approximation to the true eigenstate of -ffcffii for a finite 
value of a, this expression becomes exact in the strong- 
coupling limit, because Ji becomes very small compared 
to Jo- 



With use of Eq. (4^), En{q) is obtained as 



En{q) =ujo + Eo~ (5wo(q), 



where duj(){q) the phonon energy shift is given by 



(4.7) 



Scjoiq) -Ji—^ — V(l-cosq-a), (4.. 
A'Jo z 



with a a vector connecting nearest neighbor sites. 



Using Eqs. ( [4.lD and (|4.4|), we can rewrite Eq. ( |4.8| ) 
further as 



(5^0 (q) 



Wo 




2V, 



ph 



2U, 



ph\ 



2U; 



ph 



(4.9) 



This is our final result and it become s exac t in the strong- 
coupling limit. As explained in Sec. IV A, 5, a key quan- 
tity in this expression, depends on the nature of the 
ground state as well as the lattice structure. Specifically, 
it provides the jump in 5loo{^) at the CDW-SDW phase 
transition with the increase of U^e- Note that this ana- 
lytic result agrees quantitatively well with the numerical 
estimate for (5wo(7r) in a two-site system in Sec. ||. We 
also note that 5a;o(q) decreases monotonically with the 
increase of | Uee — 2Uph \ ■ The quantity | Uee — 2Uph IpCepre- 
sents the charge excitation gap in the atomic limiio and 
its increase suppresses the density fluctuations, leading 
to the decrease of (5wo(q). 
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V. CONCLUSION AND DISCUSSION 

In summary, we have obtained an asymptotically exact 
expression for the renormalized energy of optic phonons, 
(Do(q)j in the strong-coupling and antiadiabatic limit of 
the HH model by using the unitary transformation tech- 
niques to derive -ffeff- A jump structure in LUo{q) is found 
at the CDW-SDW regime boundary with a larger phonon 
energy shift in the SDW region. 

We have considered the insulating phase at T = 0, but 
a little different result will be obtained if T is low but 
larger than the bandwidth of the ground-state manifold, 
i.e., Jq. We treat the case in Appendix C. 

Our calculation has been done for the half-filled band. 



I f th e exchange interactions, Jg and Ji in Eqs. (3.12) and 
(4.2), are independent of the electron filling n as is the 
case for U — 0, a similar calculations can be done rather 
easily for 6luo (q) off the half-filling and we may be able 
to obtain the n dependence such as duJo{q) oc n{2 — n) 
as suggested in AC. However, in general, these exchange 
interactions depend on n, or more specifically the struc- 
ture of the state realized with these exchange interac- 
tions. This provides a highly complicated self-consistent 
problem which we have no idea to solve at present. 



where T' makes cugNph + \ U\Nd invariant, while T" vari- 
ant. In order to obtain more concrete expressions for T' 
and T", consider the polaron hopping part between i and 
j sites, J2cr ^crCjo-, in Eq- ( |3.2[ ) first. In line with Ref. |3|, 
we separate the part into three components as 



5I^tcj<T = Toii,j)+Ti{i,j)+T-i{i, 



J) 



(A2) 



Here Tq represents the processes to make Nd invariant, 
given by 

To{i,j) = ^(ni-o-cl^Cjo-nj-cr + hi^„c\^Cj„hj^„), (A3) 



with hicy = 1 — riic. On the other hand, Ti(i,j) and 
r_i(i, j) are, respectively, given by 



and 



(A4) 



(A5) 
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APPENDIX A: DERIVATION OF THE 
EFFECTIVE HAMILTONIAN 

This Appendix presents mathematical details to ob- 
tain the effective Hamiltonian H^s from H in Eq. ( ^.l[ ). 
We transform H unitarily as iJeff = e^He~^, where R 
is the antihermitian operator, determined perturbatively 
using the guiding principle that we select only such terms 
as to make LO^Nph + \U\Nd invariant among the polaron 
hopping processes in each order in t. 

Before we consider in detail, let us decompose T 
the hopping part of H into two parts as 



T = T' + T", 



(Al) 



The former (latter) describes the process which increases 
(decreases) Nd by unity. 

Now we decompose the phonon part, XjXj, in a sim- 
ilar way, but care should be taken as for the condition 
to make the energy Lo^Nph + \U\Nd invariant. Since the 
polaron hopping process is associated with the energy 
change by mU with m or ±1, wc need to separate 
XjXj into the component $fe(i,j) making the change 
in phonon number by fc = —mU/uQ and the remainder 



X/X, =$,,(», j) + j). 



(A6) 



Instead of ^kihj), we give the definition for explicitly 
as 



(A7) 



where = (a^ — aj)/V^ and J2ti\2 denotes the sum 
over £i{> 0) and ^2(> 0) except for - £2 = k. Note 
that for noninteger k — iCZ/wp, ^kihj) is nothing but 



xjXj itself. 



Combining Eqs. (|A2|) and l\A§, we obtain T' and T" 



as 



V r 



rpft 



V T" 

"1.1 



(A8) 



m=±l,0 m=±l,0 

where T/^ and are, respectively, given by 
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T'n = ~iY.^^oTm{i.m-mU/.o{i.j), (A9) 



and 



If U /lO{) is not an integer, both T[ and T^j^ vanish. 

We now perform an canonical transformation of H 
with the unitary operator as H^s — e^He^^ and 
we rewrite H^s in powers of R as 



H,s = S+[R,H] + ^[R, [R,H]] 



(All) 



where [A, B] = AB — BA. We expand i?eff and R in 
powers of t as 



H,s = J2^cS and = 
r(0) 



(A12) 



Note that W^g' = 



In first order, if we determine i?*-^) to satisfy 



if^g is simply given by 



i7« = r. 



The solution to Eq. (A13) is written as 



(A13) 



(AM) 



(A15) 



with Rm\ expressed as 



2J 



where the operator ipmihj) is defined through the equa- 
tion of 



mU cj)jn{h j) + [a\ fli + a] flj , 0m (» , j )] 



— mU/ujQ 



(A17) 



We find an explicit solution to this equation in the fol- 
lowing way: 



(atO^Ha..)'^- (A18) 



mC/-Hwo(^i -^2) 
The obtained i?^^-' is antihermitian, i.e., i?*^^^ = — i?*^^^^ 



In second order, using Eq. (A13), we can write 
H^l^ = [i?(i),r'] + [r'^^\T"]/2+[R^^\v]. (A19) 



By determining R^"^^ so as to eliminate the terms associ- 
ated with the change in ujoNph + \U\Nd, namely, choosing 
i?^^) satisfying 



(2) 

we can reduce H^^ into the following form: 

F(2)^[i?(i),r"]72, 



(A20) 



(A21) 



where the symbols prime and double prime on the square 
brackets in Eqs. (A20) and (A21) indicate that we should 



include only the processes making the energy uj^Nph 
\U\Nd invariant and variant, respectively. 



From these compact expressions for iJ^^g ±±^g , 



and h'^ , the 

forms more convenient for actual calculations can be de- 
rived. Espec ially i n the spin region, we obtain -/Jeff's as 
given in Sec. [II D . The process to reach these H^s's are 
straightforward, but a care should be taken in the reduc- 
tion of iJg^^ In general, to the process (a) in Fig.||, three 
terms out of the expression r"]'/2, i.e., T^']', 

{R^_^lTi'-rj_^R[^'^y, and {R[^^T'^^-T{'R'-_}ly contribute. 
However, in the case of ?7 > 0, only the second term 
{R^]:It{' - Ti'ii?[^V counts due to the fact that each 
site is occupied by at most one electron in the mani- 
folds with Nd = in the SDW regime. On the other 
hand, in the case of ?7 < 0, the ground-state mani- 
fold has the largest number of doubly-occupied sites and 
the condition of 'W^ = 0" in this case means that this 
continues to be the case for the excited-state manifolds 
under consideration. Thus instead of the second, the 
third term {r[^^T'!^-T{' r'i^I)' contributes to the process 
(a). At f7 = 0, no correlation works between electrons 

with different spins and all the three terms, [r''q\t^']' , 



(A16) {R^^lTi'^ry^R^'')', and {R\''T'1^-T^ R'^{)' , contribute 



(1). 



APPENDIX B: POLARON-LIQUID 
HAMILTONIAN 



As mentioned in Sec. |III Q , we can think of the region 
in the intermediate strength of a in which the polaron ki- 
netic term cannot be neglected in comparison with 



H 



(2) 
off 



In particular, we already know from the exact- 



diagonalization study on the small-size HH modelo as 
well as the perturbatiop-|theoretic study on the infinite- 
dimensional HH modeleil that "the metallic state" ap- 
pears at \U\ « 0, i.e., the transition region from CDW to 
SDW. The state is very interesting from various points 
of view includinffi-superconductivity in the alkali-metal- 
doped fuUerenesJla Thus in this Appendix, we give ffcff 
for the ground-state manifold pertinent to the investiga- 
tion of this polaron-liquid state. 

At ?7 = 0, the second -order part of H^s has already 
been given in Sec. HID. For the ground-state manifold 
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with Nph = 0, the exchange interaction Jq c an b e pro- 
vided in a more compact fashion than in Eq. (|3.8| ) as 



Jo 



-2a °° 



E 

e=i 



(Bl) 



As for the first-order part, we should include all the hop- 
ping terms, Tq, T{, and T!_-^, given in Eq. ([Ag|), because 
double occupancy is permitted at every site due to the 
vanishment of a charge-excitation energy. The sum of 
these three terms recovers the free polaron hopping term 
with the transfer integral te~". 

When \U\ is not exactly zero but small, i.e., at most 
of the order of t, w e should add the on-site interaction 
term V in Eq. ( |3.3D as it is to HeS as an additional per- 
turbation, because the correction to the operator ni^rii^ 
due to the unitary transformation of is at most of the 
order of in the ground-state manifold. Combining all 
these terms, we obtain i/cff as 

iJoff = -te"" ^ {cl^Cja + h.c.) + U ^ni-^riii 

+ 2Jo^(S,.S, + ip,p,). (B2) 

Note that this expression is valid even for n off the half- 
filling, because Jq is determined only through the phonon 
effect which is independent of n. We also note that 
the nearest-neighbor Coulomb interaction term with the 
strength V at most of the order of t, if any in the original 
model Hamiltonian, can also be added to this Hamilto- 
nian as it is. 

We shortly comment on the interesting case in which 
U/lvq happens to be equal to some nonzero integer, say, 
unity. In this situation, one phonon-excited and charge- 
excited states are degenerate. Thus, in constructing Hcs 
to describe the first-excited manifold, we should take 
account of all the one-polaron hopping terms, T^^ with 
m = and ±1, corresponding to the polaron hopping ac- 
companied with creation or annihilation of local phonons. 
Then we can expect a polaron-mobile excited state in the 
system of an insulating ground state. 



D(q, iuj^) 



2LJn 



1 



l + (C/ee/2)n(q,2W,) 



(CI) 



where 11 is the proper polarization function. Since the 
renormalized phonon energy (Dq is determined by the pole 
of -D^(q, ti^), is given by the self-consistent solution of 



1 - 



2at^on'^(q, Wo) 



l + (C/ee/2)nR(q,(I-o) 



(C2) 



where the retarded polarization function is denoted by 
n^. Thus our task is to evaluate 11^ in the antiadiabatic 
and strong-coupling {Uph 3> oJq) limit. 

In terms of real-space representation for the polariza- 
tion function I\ij{iu}u), n(q, iwi/) is expressed as 

n(q,*..,) = - e-^''-(^'-^^)]n,,(ic..), (C3) 



where we have used the relation of 11(0, iw^) — for 
Lo^ 7^ 0. Up to second order in only the polarization 
function between nearest-neighbor sites Ila contributes 
to Eq. (E|). Thus we obtain 



with 



n(q, iwi,) = - ^(1 - cosq • a)Tls,{iojy), 

a 

na(iw^) = -2t2T^[G(ie„)G'(ie„ -I- it 

n 

X A(ze„, ze„ -I- iuJi,)]'^, 



(C4) 



(C5) 



where the factor 2 in the righthand side accounts for the 
spin sum, = (2n+ l)7rT is the fermion Matsubara fre- 
quency with an integer n, G and A are, respectively, the 
one-electron Green's function and the three-point scalar 
vertex function in the atomic limit. 

For 0, the Ward-Takahashi identity relates A 

with G asE3 

.,. . , . . G~^iien + iuj^) - G-^{ien) , . 
A(ie„, ze„ + iWi,) = ■ . (C6) 



APPENDIX C: FINITE-TEMPERATURE 
GREEN'S-FUNCTION APPROACH 

In this Appendix, we consider the temperature range 
oi Jq ^ T <ti ujQ and apply the Gr een' s- function method 
to the original Hamiltonian, Eq. (2.1). In this temper- 
ature range, the electronic long-range order originating 
from the exchange interaction Jq is destroyed by thermal 
fluctuations. The incoherent nearest-neighbor electron- 
hopping processes can be treated perturbatively with re- 
spect to t from the atomic t = limit. 

From the diagrammatic consideration, the phonon 
Green's function in Eq. (^.2[) is rewritten as 



By substituting Eq. (C6) into Eq. (|C5[), we obtain Ha as 



n^iiuj^) = -At 
with x(icji^), defined by 



2x(«w^) - x(o) 



X{i^u) = -T^ G'(ie„)G(ie„ 



(C7) 



(C8) 



We can give an analytical expression for xi'i-^i') with use 
of an exact expression for jS/je„) derived through the 
Lang-Firsov transformation.E3 
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For the time being, we assume \U\ ^ T. Then xi^^u) 
is reduced into the following form in the strong-coupling 
limit at half filling: 



2Uph + \U\ 



i2Uph + \U\r-iiiu,r 



'2'Uph 



+ 



{2Uph? - {lOJ.Y 
2Uph ~ \U\ 



2 i2Uph~\U\r-{iu;,r' 



(C9) 



with / = l/(l + el'^l/^-^). Although it is not simply given 
by the limit of w,y in Eq. ( |C9| ) due to the presence 
of a singular term in proportion to l/T, x(0) can be well 
approximated by the value of limi^^_,o xl^i^!^); because 
the weight for the singular term vanishes exponentially 
for Uph > Wo- 

The analytical continuation of n(q, in Eq. ( |C^ ) 
together with Eqs. (C7) and ( |C^ ) provides n^(q, w). 
Because we are interested in the region oi uj ~ too <^ 
2Uph ± |C^|, n^(q, Lu) is essentially the same as n^(q, 0). 
Thus we can replace n^(q,cj) by n^(q,0) in Eq. 
and obtain 

5a.o(q)/wo -C/phn«(q,0), (CIO) 
up to second order in t. We can rewrite this equation as 

^wo(q) ^ V(T) f_^y ( 2UpH V 

UJO 4 \UphJ \\Uee-2Uph\+2UphJ 

x^(l-cosq-a), (Cll) 

a 

where t]{T) describes the T dependence, which is given 
by 



ry(T) 



ph 



(C12) 



up to first order in \U\/Uph- 

We have also considered the case of |[/| » T, starting 
with the inspection of x(*'^i') and obtained exactly the 
same expression for Sujq(q} as in Eq. ( pllD with / put 
equal to zero in Eq. ( |C12| ). 

At [/ = as well as for |;7| > T v(T ) is reduced to 
unity, independent of T. Then Eq. ( Cll ) pKCwides essen- 
tially the same result for Sujq as that of AC,E3l2j but it is 
different from that in Eq. (^!9|) by the factor of 1/2(1 + (5). 
This illustrates the subtlety as for the order of two lim- 
its, Uph 00 and r ^ 0. If we take the T — > limit 
first, the density fluctuations in the ordered ground state 
determines Sujo{ q) as given in Eq. ( |4.9| ). On the other 
hand, Eq. ( Cll| ) corresponds to the result with taking 
the Uph 00 limit first. In the latter limit, many de- 
generate states without the long-range order contribute 
to the thermal average even at T « 0. The absence of 



coherent density fluctuations in this case accounts for the 
reduction factor 1/2(1 + S) in SLUo{q). 
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FIG. 1. (a) An example of the results for excitation en- 
ergies in the energy range less than Zluq as a function of Uee 
with a fixed value of Uph- We obtain this result by using 
an exact diagonalization method to calculate the pole posi- 
tions in the retarded phonon Green's function D^{q,uj) in 
a two-site Hubbard-Holstein model at g = vr in the antia- 
diabatic (t — 0.2ajo) and strong-coupling {Uph ~ 2ujo) case. 
There are basically two branches for (Iiq. One is the intrinsic 
phonon mode (the sohd curve) and the other is the mode 
of electronic origin (the dashed curve). (These modes in- 
terchange their nature at Uee ~ 5u>o as indicated by the 
interchange of the solid curve with the dashed curve.) For 
comparison, the spin-excitation energy is also plotted by the 
thin dashed curve. For Uee ~ 2Uph, the third mode (the 
dotted curve), the combined charge and phonon mode, can- 
not be neglected, (b) Gorresponding spectral weight at each 
pole with the corresponding species of curve. The coeffi- 
cient of delta-function contribution in the spectral function, 
— ^lmD^{q,Lj) = Aq5{Lu — d;o(g)), is plotted &t q = n. From 
its definition, the quantity A-^ is dimensionless. 



FIG. 2. Schematic representation of the energy spectrum 
structure of the Hubbard-Holstein model in the antiadiabatic 
{t <g uoq) and strong-coupling {Uph 3> ujq) region. We can 
classify two situations, (a) The GDW {U < 0) or SDW 
{U > 0) case in which \U\ is much larger than t. (b) The 
charge-spin degenerate case in which U is about the same as 
t. 



FIG. 3. Second-order polaron hopping processes con- 

(2) 

tributing to H^^ . Two types of the processes between near- 
est neighbor sites i and j axe shown in (a) and (b), while the 
terms among adjacent three sites i, j, and k are represented 
in (c)-(e). 



FIG. 4. An example of the calculated results for the in- 
trinsic-phonon energy shift 5ujo{q) in units of Jo{U) in ID 
chains up to ten sites with the periodic boundary conditions. 
The open (solid) symbols represent the case of f/ < ((7 > 0). 
The solid and dashed curves show the dependence of 1 — cos q. 
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